Importance of second-order piezoelectric effects in zincblende semiconductors 
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We show that the piezoelectric effect that describes the emergence of an electric field in response 
to a crystal deformation in III-V semiconductors such as GaAs and InAs has strong contributions 
from second-order effects that have been neglected so far. We calculate the second-order piezo- 
electric tensors using density functional theory and obtain the piezoelectric field for [lll]-oriented 
In2;Gai_a; As quantum wells of realistic dimensions and concentration x. We find that the linear and 
the quadratic piezoelectric coefficients have the opposite effect on the field, and for large strains the 
quadratic terms even dominate. Thus, the piezoelectric field turns out to be a rare example of a 
physical quantity for which the first- and second-order contributions are of comparable magnitude. 

PACS numbers: 77.65.-j,77.65.Bn,77.65.Ly,78.67.De,71.15.-m 
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Since the discovery of piezoelectricity in 1880 by the 
Curie brothers ;l|, widespread efforts have been aimed 
at understanding this peculiar effect and developing its 
applications. Piezoelectric materials are in use today in 
a wide range of devices, including ultrasonic transducers 
for medical and sonar imaging and various types of mi- 
cropositioners and actuators. Since the early days, the 
effect has been understood as arising from displacement 
of the ions in response to a mechanical deformation, lead- 
ing to the appearance of charges on some of the crystal's 
surfaces ll]. Despite some early doubts, it is now well 
established ,2^ _3 , ^ |5|, I3 that this is a bulk eSect. It has 
two components; the contribution coming from the ionic 
displacements tends to be compensated by the purely 
electronic (frozen-ion) response, resulting in a subtle bal- 
ance between ionic and electronic contributions [31, jj, |5| . 

Until now, theoretical modeling of the piezoelectric ef- 
fect in bulk solids Q , quantum wells (Ref. |a and refer- 
ences therein) and more recently in quantum dots [sl, USL 
IllL I12I \lA |lj| has focused exclusively on the first-order 
piezoelectric tensor e^j, neglecting possible higher-order 
terms. That is, if P^ = J2j ^^mVj + \ J2jk B^ijkVjVk + •■•, 
where P is the polarization and 77 is the strain, previous 
work has concentrated on the linear coefficient e^j to the 
exclusion of the quadratic coefficient Bf^jk- 

Recent experimental determinations of piezoelectric 
constants (e.g., Refs. El El 113, El El S El) have 
tended to follow this approach, interpreting the mea- 
sured piezoelectric fields by assuming a linear relation- 
ship between polarization and strain (retaining e^j but 
neglecting B^jk)- Indeed, the experimental procedures 
used thus far to deduce piezoelectric constants from mea- 
sured fields have made it easy to overlook the importance 
of the second-order piezoelectric effect, because measure- 
ments were restricted to heterostructure quantum wells 
with small lattice mismatch (those with large lattice mis- 
match, like high-In-contcnt (In,Ga)As/GaAs structures, 
were avoided because they tend to develop unwanted dis- 
locations). Similarly, experiments for quantum wells or 



bulk materials under pressure have tended to probe only 
a very small region of strain, so that a clear signature of 
the quadratic dependence of field upon strain is difficult 
to detect. 

In this Letter we show, using self-consistent density- 
functional theory (DFT) calculations for GaAs and InAs, 
that the hitherto neglected second-order piezoelectric 
tensor gives significant contributions to the piezoelectric 
field. We show that neglecting the second-order piezo- 
electric tensor leads to an overestimation by 200% in 
the piezoelectric field for Ina;Gai_a;As quantum wells on 
GaAs in the experimentally accessible range of concentra- 
tion X. For higher In concentrations, accessible in quan- 
tum dots, second-order terms will dominate over first- 
order terms. This new insight is important because it 
represents a paradigm shift in the interpretation of mea- 
surements of piezoelectricity in quantum wells and quan- 
tum dots. 

We find it most convenient to formulate the piezoelec- 
tric response in terms of the reduced (and dimensionless) 

polarization p^ defined implicitly via Pa = fi X^uPm'^" ' 
where e is the charge quantum, fJ is the cell volume. Pa 
is the polarization in Gartesian coordinates, and a^ is 
the a'th component of the /.t'th strain-deformed lattice 
vector. We expand this reduced polarization, retaining 
the second-order strain as 
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The reduced proper piezoelectric tensor is then 
"'' " dry, - ' , 
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where we use rjj (j=l,6) to denote strain in the Voigt 
notation. Here 6°^ is the reduced proper piezoelectric 

tensor of the unstrained material, while B^jk is a fifth 
rank tensor with Cartesian coordinates fj, and the strain 
index in Voigt notation j, k and represents the first-order 



[C/vc?) as calculated from DFT. 
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change of the reduced piezoelectric tensor with strain. 
We obtain ej^ and -B^jfe from first-principles calculations 
in a manner described next. 

First-principles calculation of linear and non-linear 
piezoelectric coefficients: Symmetry considerations for 
the zincblende crystal structure imply that the only 
nonzero elements of the piezoelectric tensor are e 
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§26 = Cgg (i.e., there is only one independent element). 
Similar considerations guarantee that there are only 24 
non-zero elements of the i?^jfe tensor, which can be re- 
duced to three independent elements B114, B124, and 
-Bi56- (The other nonzero elements are obtained by ap- 
plying cyclic permutations x —> y —^ z or interchanges 
such as a: <-> y on the Cartesian axes, e.g., J5ii4 = i?225, 
-§124 = B235 = -B215, and Bi56 = i?225-) We carried 
out ab-initio calculations of these tensor elements us- 
ing a plane-wave pseudopotential approach to DFT in 
the local-density approximation (LDA) as implemented 
in the ABINIT code package|22J. First, we relaxed 
the lattice parameters for both GaAs and InAs. Next, 
linear-response calculations of the linear bulk piezoelec- 
tricconstant 65*4 were carried out on these relaxed struc- 
tures using the ANADDB module of the ABINIT pack- 
age [23, I23, which implements a direct calculation of 
the strain derivatives of the quantities of interest (Kohn- 
Sham wavefunctions, polarizations, etc.) via the chain 
rule. Then, a finite-difference technique was used in or- 
der to obtain the non-linear bulk piezoelectric tensors 
Bfijk- Specifically, we considered strain states of the form 
'7i = ?72 = % = and rj^ = rj^ ^ r]e = j for several small 
values of 7. With the strain frozen in for a particular 
value of 7, the ions were allowed to relax, after which the 
reduced piezoelectric tensor elements 611(7), ^12(7) and 
615(7) were computed using linear-response techniques 
as before. The dependence of these elements on 7 was 
then fitted, and the linear dependence extracted. From 
Eq. |(2l, this determines the three independent elements 
of the B tensor as -B114 = den/d-y, -B124 = ^612/^7, and 
-Si56 = dei^/dj. The results for GaAs and InAs are given 
in Table n 

Calculation of the piezoelectric field for large struc- 
tures (non self- consistent Poisson approach): The DFT 
method cannot be applied to lO'^-lO^-atom structures 
which are often of interest in nanoscience. Instead, such 
structures can only be calculated by non self-consistent 
methods (e.g., tight binding, k.p, empirical pseudopo- 
tentials), in which case piezoelectricity must be added 



as an external potential. Thus, to model such struc- 
tures, we first calculate the strain tensor 77 at each atom 
site using the valence force field (VFF) method. In this 
method, the bond-stretching, bond-bending, and mixed 
bending-stretching terms are derived from experiment. 
The method has been shown to give accurate atomic 
positions for defect-free bulk and alloys. For example, 
Bernard and Zunger compared strain values obtained by 
LDA and by VFF for the extreme case of a single mono- 
layer InAs superlattice in GaAs and obtaining agreement 
within 0.4% ,2Jj. From a knowledge of the strain field 
?7j(r), we can use Eq. ^ to obtain p^{y), and the piezo- 
electric charge density (per unit undeformed volume) is 
then calculated from the divergence (in undeformed co- 
ordinates) of Pfj, via 
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The calculation of the divergence [Eq. ^] is performed 
using a piecewise polynomial function to represent the 
polarization data points. Finally, the piezoelectric po- 
tential Fpiozo is obtained from a finite-grid solver of the 
Poisson equation 



Ppiczo(r) = eoV • {es(r)VFpiozo(r)} 



(4) 



where we assume an isotropic local static dielectric con- 
stant €s{y). However, note that the local polariza- 
tion cannot be defined on an arbitrarily small region of 
space 25], but only on a scale that exceeds the localiza- 
tion of the maximally localized Wannier functions [26j. 
For GaAs and InAs we average the strain tensor over 
eight-atom clusters. The piezoelectric tensor of Eq. (0) 
is position dependent since it depends on the inhomoge- 
neous strain. Furthermore, for alloys or heterostructures, 
e^j(r) also depends on the material concentration at r. 
In our case we have regions in the cell with InAs, GaAs, 
or mixed (In,Ga)As, and we use a linear interpolation 



A{y) = a;v4i„As(r) + (1 - x)Ag^a^{y) 



(5) 



of the tensors {A = e^j or -B^jfe) between the constituent 
bulk materials for the given local concentration x of the 
eight-atom cell. Finally, the solution of the Poisson equa- 
tion [Eq. Q] is obtained on the eight-atom-cluster grid 
through a conjugate-gradient algorithm with a position- 
dependent dielectric constant es(r) calculated according 
to Eq. lO with A = eg- The approach developed for the 
calculation of Vpiozo can be applied easily to very large 
(10^-atom) nanostructures. 

Testing the non-self-consistent Poisson approach: We 
tested the non-self-consistent procedure described in 
Eqs. Q-Q by comparing the results with direct self- 
consistent DFT calculations for a model quantum well 
of artificially small dimensions which can be handled by 
DFT. The system is a 30-atom [111] InAs quantum- well 



[(InAs)6(GaAs)9] epitaxially strained to the GaAs in- 
plane lattice constant. For the DFT calculations, we used 
the same pseudopotentials and convergence parameters 
as in the calculation of the piezoelectric tensors. To ob- 
tain the piezoelectric field directly from the DFT-LDA 
quantum- well calculations, we averaged out the atomic 
oscillations from the self-consistent Kohn-Sham potential 
(including ionic, Hartree and LDA exchange-correlation 
contributions). The direct DFT results of this averaging 
procedureH^ are denoted in Fig. [H as the "SCF DFT- 
LDA" curve. (A linear regression of the obtained curve 
in the region marked with "Linear Fit" gives a value for 
the piezoelectric field of 1255 kV/cm.) The piezoelec- 
tric potential obtained with the procedure of Eqs. 0}- 
0) is given in Fig. ^ as the dashed curve denoted "non- 
SCF Poisson". The potential jump at the interface is 
related to the band offset between materials j22j and is 
not present in the bare piezoelectric potential given by 
the dashed curve. The potentials are arbitrarily shifted 
to coincide in the InAs region. We see that the field de- 
duced from the non-SCF procedure is 1367 kV/cm, in 
very good agreement with the self-consistent DFT-LDA 
quantum-well result of 1255 kV/cm. 

Application to large quantum- well structures. Having 
established the validity of the non-self-consistent method 
for small systems, where comparison with DFT-LDA is 
possible, we now address nanostructures with sizes and 
composition typical of experimental conditions. We chose 
here the quantum well system Ina;Gai_a;As/GaAs where 
the In^jGai-i^As alloy is epitaxially grown on the GaAs 
substrate and the thickness of the well is around 10 nm. 
We use 12,000 atoms in the simulation cell to accurately 
represent the random alloy. We plot in Fig. [21 the piezo- 
electric potential (in mV) along the [111] direction for 
In concentration x=0.1 to 1.0. Figure [21 shows a linear 
potential along the growth axis, as expected from the 
fact that the piezoelectric charges are well localized at 
the interfaces. The small oscillations in the potential are 
due to random alloy fluctuations which are most promi- 
nent at low In concentrations. Interestingly, Fig.[21shows 
that the field reverses sign between x = 0.4 and 0.3, go- 
ing from a very strong positive field (i.e., negative slope 
along [111]) for In-rich wells to a weak negative field in 
the In-poor regime. 

The piezoelectric field extracted from the potential of 
Fig. [21 is shown in Fig. [21 as a function of the In concen- 
tration (circles). Figure[3Ib) shows the electric field with 
an emphasis on the experimentally relevant concentra- 
tion range from x= 0.10 to 0.25. We see that the ampli- 
tude of the field obtained with both linear and nonlinear 
terms is much smaller than the field obtained with e^j 
only. Considering the full concentration range, the field 
obtained using both e^j and S^^fe is shown to be nega- 
tive at low In concentrations and almost constant until 
it reaches 30%, where it reverses sign and becomes very 
strong. The field we obtain for a concentration range of 



16-20% In is nearly constant around 80 kV/cm. 

A direct comparison of calculated and measured elec- 
tric fields is difficult because only few experiments report 
the value of the measured field and the measurements 
and calculations are performed on different concentra- 
tions and thicknesses. Cho et aL|20j obtained a field of 
129 ± 12 kV/cm for a 17 % In well with an estimated 
8.7 nm thickness, and Sanchez et al. obtained a field of 
137 ± 6 kV/cm for a 17% In well with 10 nm thickness 
and 121 ± 5 kV/cm for a 10 nm thick well with 21% In 
concentration. Furthermore, the results are clouded by 
strong temperature effects [13, |2fl 1211 ( pvroelectricitv) 
and possible effects of In segregation [l9|. However, one 
experimental observation that does not seem to depend 
on well concentration and thickness is the fact that the 
measured piezoelectric field leads to an en value that is 
about 35% smaller than what is expected by using the 
linear coefficients alone. This has been reported on many 
occasions [la, llil lla, |2ij, |2S| and constitutes an unsolved 
puzzle. This result can be accurately compared with our 
calculations. If we calculate the piezoelectric field us- 
ing the experimental value of 614 for bulk InAs (—0.045 
C/m^) and bulk GaAs (—0.160 C/m^) and neglecting the 
second-order tensors B^jk in Eq. (^ (triangular symbols 
in Fig. [SI , which is equivalent to the experimental pro- 
cedure that leads to an overestimation of 614 by 35%, we 
find for the concentration region of ~18-21%, an overes- 
timation of the magnitude of the field by 34-52%. Our 
results therefore explain the ori gin of the experimentally 
observed deviation [la, [13, [iM [2il [2^ : a linear inter- 
polation between the InAs and GaAs values of 614 can- 
not reproduce the piezoelectric field of alloyed quantum 
wells since the field does not originate from the linear co- 
efficient alone but has significant contributions from the 
second-order piezoelectric tensors B^j]^ (neglected in the 
analysis of the experimental results). 

To emphasize the effect of the non-linear tensors -B^jTc 
further, we plot using square symbols in Fig.[31the piezo- 
electric field with e^j set to the DFT values and B^ij 
set to zero. The results show that, when taking only the 
linear tensor into account, the field is overestimated by 
about 200% in the region of low concentration [Fig. 3(b)], 
and even has the wrong sign at higher concentrations. 

In summary, we have shown that the second-order 
piezoelectric tensor, generally neglected so far in theoret- 
ical and experimental work, contributes significantly to 
the piezoelectric effect in zincblende semiconductors. We 
showed that the piezoelectric field calculated by includ- 
ing first- and second-order piezoelectric tensors obtained 
from DFT agree well with experiments, whereas neglect 
of non-linearities leads to qualitative disagreements. We 
argue that the "piezoelectric coefficients" that have been 
extracted from experimental work so far are actually ef- 
fective ones reflecting equally strong first- and second- 
order contributions. 

This work has been supported by U.S. DOE-SC-BES- 
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FIG. 1: (Color online) Piezoelectric potential calculated from 
Eqs. Q to Q (dashed line) and the Kohn-Sham potential ob- 
tained from self-consistent DFT calculations (solid line). Both 
calculations are for an (InAs)6/(GaAs)g superlattice epitaxi- 
ally strained on GaAs. 
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FIG. 2: (Color online) (a) Piezoelectric potential calculated 
from Eqs. 11I4II for an 11 nm Inj;Gai_a;As well epitaxially 
strained to GaAs for a; = 0.1 - 1.0. (b) for a; = 0.1 - 0.5. 

DMS under LAB 03-17 initiative and by NSF Grant 
DMR-0233925. 
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